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Abstract 

To apply lattice QCD in the calculation of glueball spectrum it is needed firstly to 
know associated operators acting on vacuum. We show how to find all the indepen¬ 
dent representations and operators, of group SO(3) PC at different dimension, since 
the work is not trivial. Then, we decompose these representation into irreducible 
representation of 0 PC group, which are listed in the note. At last we argue that 
/j(2220) and gr states can not be tensor glueball simultaneously. 
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This is a note when we calculate the glueball spectrum in quenched lattice QCD. We are 
focusing on the independent currents at each dimension here. 

As we know, in path quantization, an arbitrary bounded particle, especially, glueball, 
with J p< ~ , where J is spin of glueball, or total angular momentum of the composite particle, 
P and C are quantum number of parity and charge conjugate respectively, can be considered 
as a wave functional, which is an eigen-functional of Hamiltonian in QCD. With respect to 
the basis of eigenvalue of canonical fields, Aj = A^T a and therefore B\ = BfT a = ^ t-ijkFjk = 
2 e ijk(djAk — d^Aj — i[Aj,Ak\), a state wave functional \ij) >, particularly, glueball state with 
J PC , is generated^ by corresponding current acting on the vacuum functional |0 >, 

|V>>=0|0>, (1) 

where O is gauge invariant operator with J PC . In free QED, one can calculate operator O 
for each state. However, in QCD, one can not, due to the non-perturbational and nonlinear 
behavior of QCD. Operator O can be written as 

o = EE# 1 . (2) 

n>4 c" 

where is the i’th operator with J PC in dimension n, and c” is coefficient to solve. 
Generally, if we set single or only finite cf nonvanishing, operator O acting on vacuum 
functional will generate infinite eigenstates (physical states with definite energy) with the 
same J PC , or vice versa, there are infinite nonvanishing c” for each eigenstate. 

Therefore, to extract a definite glueball, one needs infinite operators 0\ '’s and calculates 
infinite coefficients c”’s. This is an impossible mission, since in practice we can only use finite 
operators 0- n ^ 7 s. However, thanks to functional analysis^ E], using these finite operators 
we can always obtain optimal approximation to the ultimate operator O for definite state. 
Although eigen-currents under optimal approximation also generate infinite states, we expect 
the projections of these currents on heavier states are small and we can neglect contributions 
of those heavier states. 

According to functional analysis, to extract a definite glueball, for example, quintet states 
2 ++ , we take such steps. Step one, one writes independent 2 ++ operators dimension by di¬ 
mension. In practice, the operators are surely truncated at some dimension, for instance, 
dimension seven or dimension eight. Step two, since we are working on lattice, the rotation 
symmetry is broken into cubic symmetry. We thus decompose this operators into different 
irreducible representation of O pc group. For instance, we decompose 2 +_l “ operator into E ++ 
and T£ + representation. Step three, we represent these operators on lattice by Wilson loops, 
the continuum limits of which are just the corresponding operators. Step four, for instance, 
for those nonequivalent E ++ irreducible representations belonging to different 2 ++ represen¬ 
tations in SO(3) PC group, we use variation principle to extract optimal approximations to 
different glueball 2 ++ , at finite lattice spacing a. The last step, of course, we extrapolate the 
results into continuum limit. 

This method has constraints. If we use three independent operators, which belong to three 
different representations of 2 ++ , we at most extract three different 2 ++ glueballs. Secondly, 
these operators should be independent, otherwise we shall meet singularity in the utilizing 
of variation principle. Therefore, to find all the independent currents at each dimension is 
important on lattice, since the complete basis of glueball states is generated by all the gauge 
independent currents, which, as shown, are made up of chromomagnetic fields, s, and 
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covariant derivatives, TVs, where subscripts stand for directions and superscripts are color 
indices. This is just the topic of this note. 

On one hand, for static glueball states one should sum operators over all spatial points 
on lattice to construct corresponding current. On the other hand, these operators should be 
local and gauge invariant products of B j and D{ at the same time. For instance, )F BfBf 

X 

generates 0 ++ (infinite) glueballs when it acts on vacuum (The summation over spatial points 
will be suppressed thereinafter). 

For this aim we first need to know the behavior of Bf and D bc under rotation. The 
rotation transformation is divided in two parts, one is the rotation in orbital space, and the 
other is the rotation in spin space. Bf is vector operator under rotation transformation in 
spin space and scalar operator under rotation transformation in orbital space. Therefore, 
operator Bf can be interpreted as spinor with unit spin. Whereas, D bc = di5 bc — iA^t a , 
where ( t b ) ac = if abc is generator in adjoint representation, is the vector operator only under 
total rotation transformation. This is because that the first term in covariant operator, di5 bc , 
is vector under orbital rotation and scalar under spin rotation, while the second term in 
covariant operator, —iAft a : has a reverse character, for it is scalar under orbital rotation and 
vector under spin rotation. Thus, the covariant derivative operator can not be interpreted as 
orbital angular momentum vector which ascends or descends unit orbital angular momentum. 
But, if we ignore the term —iAft a , it can. This is nearly true if the glueball is sufficient small 
and gauge fields are smooth enough, for one can always let A?(x) = 0 at one point, for 
instance the center of glueball mass, by performing a local gauge transformation. Attributed 
to this reason we sometimes call covariant derivative as orbital vector, although this will take 
many ambiguities. 

It seems firstly that finding independent current at each dimension is a trivial mission. 
However, due to the so-called Bianchi identity, one always meets some subtleties. We shall 
show how to treat these subtleties by an example, i.e. finding all the independent operators 
at the case of dimension six. Such treatment can be generalized into other dimension directly. 

All the operators at dimension six can be categorized into two types: type one is as 
Tr{(DiDjBk)Bi ) and type two is as Tr(BiBjB k ). Tr((DiBj)(D k Bi)) is not an independent 
type, for operators in such type can be rewritten as the form of type one up to a complete 
covariant derivative, which will be discarded at last. 

We first consider operators in type one. Taking advantage of standard basis of chromo- 
magnetic fields B± = q=(I?i ± and Bq = B%, we find the total spin of B k Bi can be 

0,1,2. Similarly, the standard basis D± = t(^i ± iD-i)/\/2 and Dq = D% makes that the 
total ’orbital angular momentum’ of operators DiDj be 0,1,2. It seems that there are 81 
operators. But, this is not true. 

First, we rewritten such operator using a technique of symmetry-antisymmetry decompo¬ 
sition, 

Tr((DiDjB k )Bi) = ^Tr{[(Dl j B k )B l + (DfjBi)B k ] + [{Dl j B k )B l - (D^B^Bk] 

+[(A“'£W + (DijBi)B k \ + m j B k )B l - (£>££,)£?*]}, (3) 

where Dfj = ~( D^Dj + DjDi ) and s = |( DiDj — DjDi). This technique is crucial. It 
reduces all the operators into the forms of symmetry, antisymmetry and admixture symmetry- 
antisymmetry. In each form all the operators at the same dimension make up of an invariant 
space of group SO(3) PC . Eq. @ tells us that operators in type one are divided into four 
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forms. But, terms in the last two square bracket are in fact operators in type two because of 
the identity [D t , Dj] = [F t] , ■]. (There is in fact only one independent nonvanishing operator 
€ijkTr(BiBjBk ) in the form Tr^D^-sBj^Bi).) It is easy to check the operator in the second 
square bracket is vanishing up to a complete covariant derivative. Thus, for operator belong¬ 
ing to type one it is enough to consider the operators as the form shown in the first quare 
bracket. Therefore, for operator of type one we only consider the case that the total spin and 
total orbital angular momentum are 0 or 2. There are total 36 operators to be considered 
now. 

The J PC of the total operators in type one are summarized in table 1, then, 


spin \ orbital 

0 

2 

0 

0 ++ ( i 5 0 ) 

2 ++ ( 1 D 2 ) 

2 

2 ++( 5 S 2 ) 

0++( b A)), 1 ++ m), 2 ++CD 2 ), 3 ++ ( b D 3 ), 4++(AD 4 ) 


Table 1: The J PC of the total operators in type one. The state symbols in bracket are 
obtained by interpreting covariant derivative as orbital angular momentum vector. Surely 
such interpretation has ambiguities. 

Now we turn to operators in type two. type two can be divided into two forms under 
symmetry-antisymmetry decomposition, Tr(Bi{BjBk — B^Bj)) and Tr(Bi(BjBk + B^Bj)). 
In the first form there is only one independent operator, eijkTr(BiBjBk) } which is 0 ++ . This 
operator can also be regarded as the one in type one, as shown before. We pay attention 
on the second form here, the PC of which is + —. For the case i ^ j ^ k, there is one 
independent operator, Tr(B\B 2 B 3 + B 1 B 3 B 2 ), which belongs to in cubic group and 

therefore belongs to 3 -1 in SO(3) PC group concerning to the total number of operators, 
ten, in form two. There are three operators at the case i = j = k, which construct to the 
representation T ^ . For the case i 7 ^ j = k, there are six operators. Therefore, there are 

total eleven independent operators in type two. As argued, one of them belongs to 0 ++ , 
and seven of them belong to 3 ++ . The remnant three operators belong to l -1 , since three 

operators Tr{Bi{B\B\ + B 2 B 2 + B^B^))^ make up of the representation l " 1 of the SO(3) PC 
group. 

We then have operators belong to three nonequivalent representations 0 ++ and 2 ++ , one 
representation 1 ++ , 3 ++ , 4 ++ , 1^ and 3 H respectively. 

It seems the topic is finished, but the game is not so simple, due to the Bianchi identity, 
I = D\B\ + D 2 B 2 + D 3 B 3 = 0. 

Now we consider the constraint introduced by Bianchi identity. Bianchi identity is a 
scalar operator. Then, the operators containing Bianchi identity at dimension six are as 
the form Tr(BjDjl), i.e. they belong to type one and therefore, their PC are ++• Since 
BiDj can compose representation J = 0, 1 and 2 under rotation transformation, we conclude 
that operators containing Bianchi identity can compose representations 0 ++ , 1 ++ and 2 ++ 
under SO(3) PC group. We emphasize here that for general case, we need also use symmetry- 
antisymmetry decomposition to study J PC 's of the operators containing Bianchi identity, 
although at the case of dimension six or five we need not. 

Therefore, operators belonging to 1 ++ are vanishing. The three representations 0 ++ are 
not independent and there are only two independent nonequivalent representations 0 ++ . It 
is similar for 2 ++ . For two independent representations 0 ++ , we pick on states 1 Sq and 
€ijkTr(BiBjBk ), while for 2 ++ , states 1 T >2 and 5 S , 2 respectively. 
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This statement can be checked by a directly calculation. For instance, suppose a .2 is the 
operator of state 5 Dq up to a factor, 

a 2 = Tr{-B 1 {-2D 2 + D 2 + D 2 )B l -B 2 {-2D 2 + D 2 + D 2 )B 2 -B 3 {-2D 2 + D 2 + D 2 )B 3 
+3B 2 (D 2 Di + D\D 2 )Bi + 3Bi(D 3 Di + DiD 3 )B 3 + 3B 2 (D 3 D - 2 + D- 2 D 3 )B 3 }, (4) 

and ao is the operator of state 1 S'o up to a factor, 

a 0 = Tr {B 1 {D\ + D 2 + D 2 )B 1 + B 2 (D\ + D 2 + D 2 )B 2 + B 3 (D 2 + D 2 + D 2 )B 3 }, (5) 

one has a 2 + ao = 0 if he takes advantage of Bianchi identity. 

The analysis for operators at dimension six can be generalized into other dimension. For 
instance, operators at dimension 5 belong to representations 0 h and 2 *" respectively, and 
operators at dimension 4 belong to representations 0 ++ and 2 ++ respectively.Operators at 
dimension five has no nonvanishing representation 1 *" due to Bianchi identity. 

In summary, all the operators at dimension four, five and six can be categorized into three 
independent representations 0 ++ and 2 ++ , one representation 3 ++ , 4 ++ , l -1 , S ’ 1 , 0 *" and 
one representation 2 One can make use of C-G coefficients to obtain these operators. 

These representations are irreducible in SO(3) PC group, but they are reducible represen¬ 
tations in group O pc . Since on lattice we only have O pc symmetry, we shall furthermore 
reduce these representations into irreducible representations of group O pc . The results are 
very lengthy, and we only list some of them, which will be used in our lattice simulations. 
In each irreducible representation only one operator is shown. This is enough even for repre¬ 
sentations E PC , T pc and T PC , because as long as one operator is found, other operator(s) 
belonging to the same representation can easily be obtained by suitable cubic transformations. 

• Three 0 ++ representations: 

1. The first representation 0 ++ , which is of four dimension, 

a\ oc Tr(Bj + B 2 + B 2 )- (6) 

2. The second representation 0 ++ , which is of six dimension, 

ai oc Tr[B 1 (Dl + D% + Dl)B 1 + B 2 (Dl + L% + DZ)B 2 + B 3 (Dl + D% + Dl)B 3 y, (7) 

3. The third representation 0 ++ , which is of six dimension, 


ai oc Tr[B 1 (B 2 B 3 - B 3 B 2 )\. 


( 8 ) 


Three 2 ++ representations: 


1. The first representation 2 ++ , which is of four dimension, is decomposed into irre¬ 
ducible representation E ++ and T^ + ■ For representation E ++ , 


while for T p+ , 


e 1 oc Tr{B\ - B 2 ), 
t\ oc Tr(B 2 B 3 ); 


(9) 

( 10 ) 
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2. The second representation 2 ++ is of six dimension. For representation E ++ , 

e 1 oc TriB^Dl + D 2 + D'()B l - B 2 (D 2 + D\ + D 2 3 )B 2 ], (11) 

while for T 2 + , 

t\ oc Tr[B 2 (D 2 + D\ + D 2 )B 3 }- (12) 

3. The third representation 2 ++ is of six dimension. For representation E ++ , 

e 1 oc TrlBtiD 2 - D 2 )B l + B 2 (D 2 - D 2 )B 2 + B 3 (D 2 - D 2 )B 3 ], (13) 

while for T 2 + , 

t\ oc Tr\B\D 3 D 2 Bi + B 2 D 3 D 2 B 2 + B 3 D 3 D 2 B 3 ], (14) 

• Representation 4 ++ only occurs in dimension six. This representation is reduced into 
A^ + , E ++ and two T-^+’s. For representation Af + , 

a\ oc Tr[B l (2D 2 l -Dl-Dl)B 1 + B 2 (-D 2 l +2Dl-Dl)B 2 
+B 3 (—D\ — D\ + 2D 3 )B 3 — 2B 2 (D 2 Di + DiD 2 )B\ 

—2B 1 (D 3 Di + DiD 3 )B 3 — 2B 2 (D 3 D 2 + D 2 D 3 )B 3 ], (15) 

and for representation E ++ 

e 1 oc TrlB^Df - Dl)B 1 - B 2 (D 2 - D 3 )B 2 - B 3 (D\ - D 2 )B 3 

—4BiD 3 D\B 3 + 4B 2 D 3 D 2 B 3 ], (16) 

• Representation 3 ++ , which is reduced into A 2 + , T 1 ++ and T 2 + , only occurs in dimen¬ 
sion six. For representation A 2 + , 

a 2 oc TrlB^D 2 - D\)B X + B 2 (D 2 - D 2 )B 2 + B 3 (D 2 - D 2 )B 3 ], (17) 

and for representation 

t\ oc Tr[—B 2 D 2 D 3 B 2 + B 3 D 2 D 3 B 3 + 2B 2 D\D 3 B\ 

—2Bi(D 2 Di + D\D 2 )B 3 + B 2 (D 2 — D 2 )B 3 ], (18) 


We meet representation 0 + in dimension 5, 

ai oc Tr[B 2 D\B 3 + B 3 D 2 B\ + B\D 3 B 2 ]. 

Representation 2 h is met in dimension 5. For E 

e 1 oc Tr[B 2 D\B 3 + BiD 2 B 3 ], 

and for T 9 _+ 


(19) 

( 20 ) 
( 21 ) 


t\ oc Tr[B\D 3 B 3 + B 2 D 2 B\\. 

Representation l -1 at dimension 6 is reduced into irreducible representation T ^ , 

t\ oc Tr[B\B\ + B\B 2 + B\B 3 ]. (22) 
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( 23 ) 


• Representation 3 + occurs in dimension six. We show the operator in A£ , 


a 2 oc Tr[B 2 BiB 3 + B 2 -B 3 R 1 ]. 


Now we have listed the results of decomposing twelve independent representations into 
irreducible representations in O pc group. If we use variation principle we obtain twelve 
glueballs. these glueballs, except excited scalar and tensor glueballs, have been measured by 
other authors. The listed results tell us the structures of different glueballs, which is useful, 
for instance, in glueball decay or in condensate of different current in glueballs. 

Although we do not make any lattice simulation here, we can obtain some physical con¬ 
clusions from the results, one of which will be shown in the next paragraph. 

to apply variation principle we can obtain 12 glueballs, which are in fact the optimal 
approximation of the real glueballs. Attributed to dimension analysis, it is reasonable to 
regard 4 ++ glueball is the heaviest glueball with PC = ++ among those 12 glueballs. As 
pointed out by previous works®. the mass of which is about m 4 ++ = 3 ~ 4 Gev. We have 
shown elsewhere that it is enough to use chromomagnetic fields and their covariant derivative 
in the calculation of glueball spectrum merely. Therefore, there are three and only three 
2 ++ states lower than m 4 ++. One possibly debates that their masses will be corrected by 
correction from fermion-loop and from higher-dimension operators. On one hand, to consider 
the correction of fermion-loops we should go into full QCD. However, it is believed that the 
quenched approximation almost does not change the ratio between different states. We think 
the mass correction from fermion-loop is mainly on scale shift. On the other hand, we shall 
get decreased glueball masses on lattice if we consider the correction from operators at higher 
dimension. Such decrease is both on 4 ++ glueball and on 2 ++ glueballs. We conclude that 
mainly effect of higher dimension operator is also a scale shift. Therefore, the statement, there 
are three and only three 2 ++ states lower than m 4 ++, is right or almost right even including 
the two corrections, possibly up to a scale shift. If fj( 2220) (or £(2230)) is a tensor glueball, 
one must find other two tensor glueballs with mass lower than m 4 ++. But, since there are 
only three tensor glueballs below m 4 ++, it is difficult to place three gx states, which are more 
like tensor glueballs. However, if gx are interpreted as three glueballs, fj{ 2220) should never 
be tensor glueball. The tensor glueball interpretations of fj (2220) and the tensor glueball 
interpretation of gx are conflictive interpretations, therefore. Since there are more supports 
of the tensor glueball interpretation of gx, at least in nowadays, fj( 2220 ) can not be tensor 
glueball, if it exists. 

There are quenched lattice gauge calculation supporting the tensor glueball interpretation 
of fj(2220 ), for instance, reference f5J claims that masses of ground scalar and ground tensor 
glueball are 1.73 Gev and 2A0Gev respectively. The scale parameter is determined by phe¬ 
nomenal Sommer potential 0 in these works ®, that is , ro is determined by r 2 (^| r= o. 49 / m = 
1.65, since there is no exclusive experiment result on glueball mass. However, ro dependents 
on different phenomenal potential models, for instance, as author pointed out in reference | 5 ], 
the logarithmic potential gives values that are lower by about 10%. Combining the argument 
in the last paragraph, we think that there may exist a global scale shift in these works ^ ^ . If 
this is the case and the ground tensor glueball is gx( 2050), the mass of scalar glueball should 
be 1.73 x x^Gev = 1.48 Gev, which looks very like /o(1500). If the ground scalar glueball is 
/o(1500), the dominant component should be ss in /o(1710). This interpretation is favored 
in many phenomenon analysis. For instance, we obtain a more natural conclusion that the 
ss — nn mass difference is about 300 MeV ^ . In fact, a more recent work® reveals that one 
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can not distinguish the scalar glueball interpretation of /o(1500) and /o(1710) nowadays in 
unquenched lattice QCD. From the above argument we also obtain the mass of 1 glueball 
should be 3.85 x § ^Gev = 3.29 Gev. This means that 1 glueball mass and J/ip mass are 
very close. Therefore, taking account into the lattice computing error, to solve the ” p — ir 
puzzle” of «//?/> and ip 1 decays, the suggestion^ ® that J /ip —► pir is enhanced by a mixing 
of the J/ip with an 1 glueball that decays to pir is also suitable. 

Unfortunately, in simulations we found that the mass of first excited tensor glueball and 
the mass of second excited tensor tensor glueball are very close. Such unexpected occasion 
makes the traditional approach of variation principle, illustrated by C. Morningstar, has 
disadvantages in the calculation. We will show an alternative method in the future. 

In summary, we have finished the first two steps in the calculation of glueball spectrum. 
The leavings will be shown in the future. Although we do not make any lattice simulation 
here, we obtain an interesting conclusion, that is, fj(2220 ) and gx states can not be tensor 
glueball simultaneously. 
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